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^ , Abstract 

^ ■ I suggest wave equations for the scalar, pseudoscalar, vector, and 

. pseudovector fields with different masses for spin zero and one states. 

\ Tensor, matrix, and quaternion formulations of fields with two mass 

' and spin states are considered. This is the generalization of the Dirac- 

- Kahler equation on the case of different masses of fields with spin one 

^ I and zero. The equation matrices obtained are simple linear combina- 
tions of matrix elements in the 16-dimensional space. Spin projection 
operators and solutions of equations (for spin one) in the form of 

, matrix-dyads are obtained. The canonical quantization of fields un- 

I der consideration is studied. The anomalous interaction of the scalar, 

2 ■ pseudoscalar, vector, and pseudovector fields with the external elec- 

, tromagnetic field is considered. Three constants which characterize 

I the anomalous magnetic moment and quadrupole electric moment of 

Q^' a particle are introduced. 
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1 Introduction 



^ I The modem high energy physics require the consideration of fields with 

higher spins. So, supersymmetrical (SUSY) field models introduce super- 
partners of ordinary particles which are high spin particles. Cosmological 
models also deal with fields possessing different spins. 

The theory of relativistic wave equations describing particles with arbi- 
trary spins has a long history p, but is not complete yet. The simple high 
spin fields are fields with spin one. Ordinarily, spin-1 particles are described 
by Proca j2] or Petiau-Duffin-Kemmer equations [2], [HI- Anomalous in- 
teractions of vector particles with electromagnetic fields were considered in 
[Oj , [Z| , [H] • Other schemes to the theory of spin one particles were developed 

by i, m- 
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Kahler jTTj considered an equation for inhomogeneous differential forms 
which is equivalent to a system of scalar, vector, antisymmetric tensor, pseu- 
dovector and pseudoscalar fields [12], [IH]. Now, the Kahler equation (called 
Dirac-Kahler equation), in the framework of differential forms is intensively 
used for description of quarks on the lattice Note that some authors J3] 
considered equations for the system of antisymmetric tensor fields which are 
equivalent to the Kahler equation before the appearance of [TJI. It should 
be mentioned that the Lagrangian of a system of scalar, vector, antisymmet- 
ric tensor, pseudovector and pseudoscalar fields (equivalent to Dirac-Kahler 
fields) is invariant under the internal symmetry group 5*0(4,2) (or locally 
isomorphic group SU{2,2)) [12], [I3]- 

In this paper, we generalize the Dirac-Kahler equation on the case that 
scalar, pseudoscalar fields having different masses, mo, compared with vector, 
pseudovector fields, m. At the particular case mg = m, one arrives at the 
Dirac-Kahler fields. 

The paper is organized as follows: In Section 2, we introduce the field 
equations which are the generalization of the Dirac-Kahler equations. It is 
prove that fields introduced describe scalar, pseudoscalar fields with mass 
mo and vector, pseudovector fields with mass m. We show that introduced 
16-component system of equations can be represented as two 8-component in- 
dependent sub-systems with self-dual, and antiself-dual antisymmetrical ten- 
sors. At parity transformations, these systems are converted to each other. 
In Section 3, we construct the matrix form of 16-component relativistic wave 
equation. It is shown that matrices of the equation obey the Dirac commuta- 
tion relation. This equation also contains two projection operators which are 
connected with two possible mass states. Section 4 is devoted to the matrix 
form of 8-component equations. In the case of equal field masses, rriQ = m, 
the 8-component matrices also satisfy the Dirac commutation relation. In 
section 5, we study the Lorenz covariance of equations under consideration. 
The one-parameter symmetry group of the Lagrangian and the correspond- 
ing conserve current are investigated here. The quaternion form of fields is 
constructed in section 6. This quaternion form includes different particu- 
lar cases of field schemes such as Proca, Maxwell, etc. In section 7, spin 
projection operators and solutions of equations for spin one in the form of 
matrix-dyads are found. The quantization of fields is carried out in Section 
8. In section 9, we consider the possible anomalous electromagnetic field 
interaction. Section 10 contains the discussion. The appendix is devoted to 
the computation of products of 16-component matrices. 
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The system of units h = c = 1, a = e^/in = 1/137, e > is used. 

2 Field equations 

Now, I suggest new field equations for the set of bosonic fields which have two 
spin one states and two spin zero states with differen masses. In particular 
cases, these equations are transformed to one considered in literature. 

Let us introduce and investigate the system of equations for free scalar 
ip{x), pseudoscalar <^(x), vector V5^(x), pseudovector <^^(x) and antisymmet- 
ric tensor f^ui^) fields: 

d^Lp^^{x)-d^Lp{x)+m'^Lp^{x) = 0, d^(f^^{x)-df,(p{x) + m'^(ff,{x) = 0, (1) 

ad^Lpf,{x) - (^{x) = 0, ad^(ff,{x) - <^(x) = 0, (2) 
<^fiu{x) = d^ipu{x) - dyLp^{x) - e^uapda'f'iiix), (3) 
where = {dm, —ido), x^ = {xm,ixo)] Xq is the time, and 

'?fiu{x) = ]^ef,yaP'^al3{x) (4) 

is the dual tensor, e^yap is an antisymmetric tensor Levy-Civita; £1234 = —i] 
a and m are parameters which are connected with the masses of scalar, 
pseudoscalar and vector, pseudovector fields. Equations (l)-(4) are the gen- 
eralization of the tensor form of the Dirac-Kahler equation (TJ. At a = 1, 
we come to Dirac-Kahler equations [12] (see also [H]); the case ^{x) = 
and if^ix) = was considered in [1^]; at a = 1, (p{x) = 0, <^^(x) = 0, we 
arrive at the Stueckelberg equations jU] (see also ^7\); the case a = 0, cor- 
responding vector and pseudovector fields, was investigated in TF; in the 
case ip{x) = 0, (^(x) = 0, <^^(x) = one arrives at the Proca equations [2]; 
at a = 1, ip{x) = 0, (pfj.{x) = 0, m = 0, we obtain generalized Maxwell's 
equations investigated in [1^]. So, field equations (l)-(4) include many im- 
portant cases of equations for bosonic fields with different number of freedom 
degrees. 

To clear up the physical meaning of constants a, m, one can get from 
Eqs. (l)-(4) the second order equations for vector and pseudovector fields 

dlip^{x) - (1 - a)df,du^u{x) - m^ip^{x) = 0, 

(5) 
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dl(p^{x) - (1 - a)d^dyip^{x) - m^ip^{x) = 0. 

Because two equations in (5) are identical, masses of vector and pseudovector 
fields are the same. To obtain the masses of vector (and pseudovector) fields, 
we write the equation for vector field (5) in momentum space (see [TH]): 

(pI + m^) 6^^ - (1 - a)pf,p^^ Lp^{p) = 0, (6) 

where p"^ = pj, = p"^ + pi = p"^ — Po- The matrix 

M={p^ + m^) J4 - (1 - a) (p.p) , (7) 

where the I4 is the unit 4x4- matrix and the matrix-dyad {p.p) has the matrix 
elements (p-p)^^ = P^Pv, obeys the equation 

[m -rn^ - [m - - ap^] = 0. (8) 

It follows from Eq. (8) that the matrix M possesses the eigenvalues Ai = 
p^ + m? , A2 = ap^ + rn^ ■ So, Eq. (6) has nontrivial solutions at detM = 
or Ai = A2 = 0. This leads to the mass spectrum of the vector field ^u{.x): 

2 

on p TTlj 

p = —m , p = . (9) 

We can draw a conclusion from Eq. (9) that one state of the field ipu{x) 
corresponds to the mass m, and the second state corresponds to the squared 
mass 

ml = — . (10) 

(Jj 

The requirement that the masses are real (m^ > 0, ml > 0), we get the 
condition a > 0. To identify the the mass states with the spin states, one finds 
from Eqs. (l)-(4) the second order equation for the scalar and pseudoscalar 
fields: 



m 



2 



dtifix) v{x) = 0, 



« (11) 

d'ivi.^) - — <^(a;) = 0. 
a 

From Eqs. (10), (11) we conclude that scalar '{){x) and pseudoscalar (p{x) 
fields with spin s = possess the mass m^. Therefore, the vector and pseu- 
dovector states with spin s = 1 have the mass m. In the case of a = 1, states 
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with spin s = 1 and s = possess the same mass m. For arbitrary variable 
a > the system of equations (l)-(4) describes fields which may have the 
vector and pseudovector states with the mass m and scalar and pseudoscalar 
states with the mass mo- The Dirac-Kahler equation [llj is obtained by 
setting m = rriQ (a = 1). 
Introducing the variables 

M{x) = {<^{x) - i^{x)) , N{x) = ((^(x) + i^{x)) , 

= ^ i^f^i^) + i^t^i^)) , Nf,^{x) = i^^,Ax) + i(p^,^{x)) , 
and adding and subtracting Eqs. (l)-(4), we obtain equations 

dMi^Ax) - d^M{x) + m^M^{x) = 0, a(9^M^(x) = M{x), 

M^^{x) = d^M^{x) - dyM^{x) - iei^^^pdaMpix), 
duN^uix) - df,N{x) + m^Nf,{x) = 0, adf,N^{x) = N{x), 

N^u{x) = d^N^{x) - d^N^,{x) + ie^yapdo,Np{x). 



(12) 



(13) 



The tensor M^y(x) is self-dual tensor, M^^{x) = —iM^^{x), and the tensor 
N^^{x) is the antiself-dual tensor, N^y{x) = iN^^^x). 

The self-dual tensor M^j, possesses 3 independent components and real- 
izes the (1, 0)-representation of the Lorentz group. The antiself-dual tensor 
N^u also has 3 independent components and is transformed under the (0, 1)- 
representation of the Lorentz group. Eqs. (12) and (13)) are 8-component 
equations which describes eight independent variables (M(x), My{x) , Mah{x)), 
{N{x), N^{x), Nab{x)) are not invariant under the parity transformations 
separately. At P-transformations, Eqs. (12) are transferred into Eqs. (13) 
and vice versa. It is possible to have Lagrangian formulation only for the 
whole system of equations (12), (13) or (l)-(4). P-noninvariant equations 
(12) (or (13)) have no the Lagrangian formulation. Eqs. (l)-(4) realize the 
(0, 0) © (1/2, 1/2) © (1, 0) © (0, 1) © (1/2, 1/2) © (0, 0)-representation of the 
Lorentz group, and are equivalent to equations (12), (13)). 
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3 Matrix form of 16-component equations 

Now we obtain the matrix form of equations (l)-(4). To get the relativistic 
wave equation in the matrix form, we introduce new variables 

^'/.(x) = irrv^f,{x), tpo^x) = -i<^{x), e^ua/s = ie^^af} (ei234 = !)• 
Then Eqs. (l)-(4) with the help of Eq. (10) are rewritten as 

2 

df,i^^,{x) + ^M^) = 0, 



m 



^u'ip[^,u]{x) + d^ipo{x) + mtpf,{x) = 0, 
dui^fi{x) - df,i)u{x) - e^uapdai'isix) + rmp[^i,]{x) 



df^ij^ix) + -^M^) = 0. 

m 

With mo = m, we arrive at the Dirac-Kahler equations. It is convenient to 
introduce the 16-component wave function 



(14) 



^'(x) = {ipAix)} 



A = 0,/i, [fxu],fx,0, 



(15) 



where ip- = ip^, ip^ = ipo- Let us introduce the the elements of the entire 
algebra e^'^ 20^ with dimension 16 x 16. The product of these matrices and 
their matrix elements are given by 



CD 



^AC^BD, 



(16) 



where A,B,C,D = 1,2,..., 16. With the help of the matrices e"^'^ and Eq. 
(15), equations (14) become 



+ 7:e 



+ 



+ e 

2 



[puj] ,p 



+ 



AB 



(17) 



2 J m 



J AB 



We note that the matrices 

are the projection matrices which obey the equations 

PoPl = PlPo = 0, Pi+Po = /l6, 

where the Jig is the unit 16 x 16— matrix. Introducing the matrix 

r, = e^'I^'^l + el'^^]'^ + e'^'O + e°'^ + e^^ + e^^ + ^e^,^^ (e^'IH + ^M,?) ^ (19) 
Eq. (16) becomes 

V,d, + mPi + ^Po ^(x) = 0. (20) 
m J 

Eq. (20) is the matrix relativistic wave equation which describes the system 
of the vector (pseudovector) and scalar (pseudoscalar) fields with the masses 
m and mo, respectively. The 16 x 16— matrix Ty can be cast as follows: 

^(1) _ + ^[mH,/^^ ^(T) = Ig^^^^ + ^ (21) 

The 10- dimensional matrices /S^^-*, (3^^ and 5-dimensional matrices /J^^-*, (3^'^ 
obey the algebra [3], jl], 

PfxPuPa + PaPuPf, = S^uPa + SauPfi- (22) 

The 16- dimensional matrices Z?^"-* realize the reducible representations 
of the Petiau-Duffin-Kemmer algebra (22) [21]. The 16 x 16-matrix T,, sat- 
isfies the Dirac algebra: 

r,r^ + r^r, = 25^,. (23) 

In the case when the masses of all states are the same, mo = m, Eq. (23) is 
transformed into the matrix form of the the Dirac-Kahler equations ^2j : 



(r,9, + m)^(x) = 0. (24) 
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4 Matrix form of 8-component equations 

To get the matrix form of 8-component equations (12), we introduce four- 
component wave functions: 

where Ma{x) — {l/2)€amnMmn{x). With the help of Eqs. (25), the system of 
P-noninvariant wave equations (12), may be cast in the form 

(26) 

where we introduce matrices: 
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(27) 



(28) 



It is easy to verify that four-dimensional matrices (a = 1, 2, 3) obeys the 
Pauli commutation relations: 

{ai,ak} = 25ik, afc] = ^ieiki^i, (29) 
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where €123 = 1. Note, that matrices may be obtained from matrices 
by setting a = 1, and = —1!^. Equations (26) can be also written in the 
form of one 8-component equation: 



f3^d^ip{x) + mip{x) = 0, 



(30) 



where the matrices of the equation and 8-component wave function (p{x) 
are given by 



ax) 









(31) 



In the case a = 1 when the masses scalar and vector states equal each other, 
m = mo, matrices (at a = 1) obey the Dirac commutation relations (23) 

(see na, m)- 

In the same manner, to obtain the matrix form of 8-component equations 
(13), we introduce four-component wave functions: 



ax) 



-im 



Na{x) 

Ni{x) 



X {X) 



Na{x) 

N{x) 



(32) 



where Na{x) = {1 / 2)eam.nNmn{x) . Then Eqs. (13) can be represented by 

Oi'^,df,aix) = mx'ix), 



(33) 



^'Ax\x) = mCix), 



where 



a 



f i \ 

-i 

i 

V m y 

/ -i 0\ 

z 

i 

V ia J 



an 



a 



/ i 0\ 

i 

-i 

V ia J 

f -i \ 

-i 

-z 

V ia J 



(34) 
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0l^ 



/ -i\ 



i 

V -z / 

/ -i \ 

i 

-i 

V -i y 



/ z 

-z 


V -i / 

^ i 

i 

i 

V / 



Two equations (33) are rewritten by one 8-component matrix equation 

/?;9^</.'(x)+my.'(x) =0, (35) 
where matrices and columns '^'[x) are 




a, 



(36) 



The matrices f3'^ also at a = 1 obey the Dirac algebra (23) (see |22]- One 
may combine Eqs. (30), (35) in the 16-component relativistic wave equation, 
as follows 



(n^9^ + m)^(a;) = 0, 



where 



i^{x) 



n, 




(37) 



(38) 



The 16x16 - matrices 11^ obey at a = 1 the Dirac algebra (23). The whole 
system of equations (12), (13) is equivalent to Eq. (37) or Eq. (20). Equa- 
tion (12) (and (13)) as well as Eq. (30) (and (35)) are parity noninvariant 
separately and at the same time the system of the two equations (12), (13) 
(or Eq. (37)) are P-invariant. 

We note that at the particular case m = niQ (a = 1), Eq. (30) and 
Eq. (35) are invariant under the GL{2,c) symmetry group [22]; [IHI- In 
this case Eq. (20) and Eq. (37) are equivalent to the Dirac-Kahler equation 
and possess the G'L(4, c) internal symmetry group, but the Lagrangian is 
invariant under the 50(4, 2) group [12], [IH]- The condition m ^ mo (a ^ 1) 
eliminates the degeneracy and destroys the symmetry. 
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5 The Lorentz covariance and symmetry of 
16-component wave equation 

At the Lorentz transformations of the coordinates: 

— L^iiX^^ (39) 
the matrix L — {L^j,} obeys the equation 

LfiaLua — ^fiu- (40) 

The wave function (15) under the Lorentz coordinates transformations (39) 
is converted into 

^'(x') =T^(,t), (41) 

where the T is the 16-dimensional reducible tensor representation of the 
Lorentz group. Taking into account that at the Lorentz transformations (39) 
the derivatives 9^ transforms as 9^ = L^^dy, the first order wave equation 
(20) becomes 



r,a:, + mPi + ^PoU'(a;') 



(42) 



r^L^yd^ + mPi + ^Po j T-^ix) = 0. 
Eq. (20) is form-invariant if the equahties 

r^TV^rr,, p^t^tp,, p^t^tp^ (43) 

are satisfied. Let us use the infinitesimal Lorentz transformations (39) with 
the matrix 

^jxv ^jxv ~l~ ^jxvi ^ jxv ^v^xi (44) 

and the infinitesimal transformations (41) with the matrix 

T = 1 + i^^. J^., (45) 

where 1 = /le, J^v are the generators of the Lorentz group in 16— dimensional 

space, are six parameters defining rotations and boosts. Replacing equa- 
tions (44), (45) into Eqs. (43) and taking into consideration the smallness of 
parameters we arrive at the equations 
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It is easy to verify with the help of formulas of Appendix A, that generators 

J,u = \{r^,ru] + ri,T,]), (47) 

satisfy Eqs. (46). The matrices Ti, are given by [21] 

r, = /3W-/5(-), (48) 

and also obey the Dirac algebra (23) and commute with the matrices F^: 



r r = r r 



Generators of the Lorentz representation (47) are also valid in the case of 
m = mo 5 QHI- The 16-dimensional matrix of the finite transformations 
(41) is given by 

T = exp (le^,J^,) . (49) 



.2 

Let us consider a relativistically invariant bilinear form 

'*(x)^'(x) = ^+(x)r7^f(a;), (50) 

where is the Hermitian-conjugate wave function and t] is the Hermitianiz- 
ing matrix in 16-dimensional space. As for the case of Dirac-Kahler equation 
in the matrix form |2I], ^2], ^Bj, it is given by: 

V = r4r4. (51) 

The matrix r] obeys the necessary (see, for example j^O], [231) equations 

r]Ti = -Til] (i = 1, 2, 3), r]Ti = T^r]. 

Now we consider the internal symmetry group of Eq. (20). With the help 
of Eqs. (16), (18), (21), (48) it is easy to show that the matrix = rir2r3r4 
commutes with the operator of Eq. (20): 

Y^d^ + mPi + -^Po j = f r^5^ + mPi + -^Po j Y^. (52) 

This means that the matrix Fs can be considered as the generator of the 
symmetry group GL{\, c) of Eq. (20): 

^'(x) =exp(cjr5)^(x), (53) 
12 



where uj is the complex group parameter. The Lagrangian corresponding to 
the first order relativistic wave equation (20) is given by 



-4 



^(x) r^9^ + mPi + ^Po ^(x) 
V m 



(54) 



The Lagrangian (54) is invariant under the transformations (53) at the re- 
striction: uj* = 00, i.e a; is a real parameter. This constrain leaves GL[1,R) 
subgroup of the Lagrangian symmetry. The transformations of the Lorentz 
group (49) commute with transformations (53) of the internal symmetry 
group. As a result, the parameter of this group is a scalar. 

According to Noether's theorem, the invariance of the Lagrangian (54) 
under the transformation (53) provides the conservation tensor (see also [T^ . 
[THj for the case m = mo): 

K^ = ^{x)T;r,^{x), d^K^ = 0. (55) 

The subgroup U{1) of gauge transformations: ^'(a;) = exp (za) \E'(x) {a* = 
a) leads to the conservation of four-current = i'^{x)Tf^'^{x): d^J^ = 0. In 
the case of the degeneracy m = tuq, the symmetry group of the Lagrangian 
(54) is the ^0(4,2) group [HI, 

Using the formulas of Appendix A, we obtain the matrix = rir2r3r4: 

= - - e^'^ - e>'~^ - ^e^^p^e^^"'^'^'^^ (56) 

With the aid of Eqs. (15), (56) the infinitesimal transformation (53): (5\E'(x) = 
rr,'${x)6uj leads to the transformations of the fields (15): 

SiJo{x) = -tjjo{x)6u, 6%[)q{x) = -tpo{x)6uj, Stp^i^x) = -ij^{x)Suj, 

1 

Slpf,{x) = -1p^{x)SuJ, Si![^^]{x) = --e^ugatp[ea]ix)5uj. 

The finite transformation (53) can be represented as 

\E''(x) = (coshcij + Fs sinha;) \E'(x), 



which is an analog of dual transformations in electrodynamics but with the 
imaginary value of the angle lu [TH] • 



13 



6 Quaternion form of field equations 

Now we obtain the quaternion form of equations (l)-(3). Quaternions may 
be considered as a generalization (doubling) of the complex numbers (see, 
for instance, |21])- Quaternions are useful tools for analyzing the symmetry 
of fields. 

The quaternion algebra consists of four basis elements = (efc,e4) (see, 
for example [24]) with the products: 

el = 1, el = el = el = -1, 6162 = -6261 = 63, 6263 = -6362 = Ci, 

(58) 

6361 = -6163 = 62, 646^ = 6^64 = Cm [m = 1, 2, 3), 

where 64 is the unit element. The biquaternion (or complex quaternion) q is 
defined as 

Q = Qf^ef, = Qmem + ^464, (59) 

where are complex numbers. With the help of Eqs. (58), the product of 
two quaternions, q, q', is given by: 

QQ' = Ma - Qmq'm) 64 + {qigm + QaqL + emnkqnq'k) em- (60) 

It should be noted that the combined law for three quaternions: (^1(72) = 
Qi (92^3) holds. The quaternion conjugation (hyperconjugation) is defined as 

q = ^464 - gmCm = 94 - q, (61) 

with the equalities qi + q2 = qi + ^25 ^M2 = ^2^1- 

To get the quaternion form of Eqs. (1), we rewrite them as follows: 

doEmix) + dm^{x) - emnkdnHk{x) - m^ip^{x) = 0, 

d^Em{x) + dof{x) + m^ip{x) = 0, 

(62) 

doHmix) + dmf{x) + emnkdnEk{x) - m^(pm{x) = 0, 

dmHm{x) + dQLp{x) + m^(^(a;) = 0, 

where ^^{x) = {ipm{x) , iipo{x)) , Hm{x) = ^emnk'^nk{x) , Era{x) = iiPmAix), 

do = dt {t is the time). Introducing the biquaternions 

V = Cf^df,, F = F^e^, G = G^e^, 
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(63) 

with the help of Eqs. (60), we represent Eqs. (62) as follows 

VF + m^G = 0. (64) 

So, quaternion equation (64) is equivalent to Eqs. (1). To represent Eqs. 
(2), (3) in quaternion form, one may cast them into 

ip{x) = a {doipo{x) + dm^m{x)) , ^p{x) = a {do(po{x) + dm'f>m{x)) , 

Hm{x) = -dm(po{x) - do(pm{x) + emnkdn^k{x), (65) 
Em{x) = -dm^Po{x) - doipm{x) - emnkdn^k{x). 

It easy to show with the help of Eq. (60) that equations 

VG = d^G^Ci + {diGra — draGi — emnkdnGk) Cm, V = 6^9^, 
= d^G^Ci — (194^^ — draGi — emnkdnGk) Gmi Gfi = (^4, —Cm) 

are valid. The arrow above the V shows the action of the differential operator. 
Then Eqs. (65) take the form 

F = i^G*^ - 1±^VG. (66) 

At the case a = 1 {m = mg), one arrives at the quaternion form (Eqs. (64), 
(66)) of Dirac-Kahler equation [TB] . 

Under the the Lorentz group transformations, the quaternions G, F, V 
and V are converted into (see [12], [H]) 



G"^ = L GL, = LFL, 



(67) 



= L*VL, = LVL*, 

where quaternions L obey the constrain LL = LL = 1. The field equations 
(64), (66) are invariant under the Lorentz transformations (67). 

From Eqs. (64), (66) one can obtain the quaternion form of some equa- 
tions for the vector fields (the Proca equations. Maxwell's equations and so 
on). 
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7 Spin projection operators and density ma- 
trix 

To obtain all independent solutions of the matrix equation (20) in the form 
of matrix-dyads, one needs to construct projection operators [22] extracting 
"pure" states. With the help of generators of the Lorentz group (47), we find 
the spin projection operator (see [12], [THj): 



= ~T71 \^abcPaJbc = i^^abcPa {'^b'^c + ^b^c) (68) 

2|p| 4|p| ^ ^ 

This operator obeys the "minimal" equation as follows: 

ap{cxp-l){ap + l) =0. (69) 



According to the general method one obtains the projection operators 
extracting spin projections ±1 and 0: 



S(±i) = {ap ± 1) , 5(0) = 1 - aj. (70) 

These operators satisfy the ordinary relationships: •S'^^-ti) = S(±i)S(o) = 

0, S(Q) = 3(0). We consider fields with different spins, one and zero. To 
separate these states, one has to construct the squared spin operator. Such 
an operator is given by the squared Pauli-Lubanski vector a^: 

One can verify with the aid of Eqs. (47) (see Appendix) that this operator 
satisfies the matrix equation 

a' (a^ - 2) = 0. (72) 

Eq. (72) shows that eigenvalues of the squared spin operator cr^ are s(s + l) = 
and s(s + l) = 2; that corresponds to spins s = and s = 1. The projection 
operators separating these states are given by 
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Operators (73) possess the properties 

The projection operators S'^q^, S"^-^^ acting on the wave function extract pure 
states with spin and 1, respectively. 

The system of field equations under consideration includes two fields with 
spin one (vector ipfj,, and pseudovector ipf^, fields), and two fields with spin zero 
(scalar ipQ, pseudoscalar ipo, fields). So, there is a doubling of the spin states 
of fields, i.e. degeneracy. It is necessary to introduce additional projection 
operators to separate these states and to have pure spin states. For this 
purpose, one may explore the internal symmetry of fields which mixes scalar 
and pseudoscalar fields as well as vector and pseudovector fields, and is given 
by Eq. (53). We introduce the projection operator as follows: 

Ax = 1(1 + XT,) , (74) 



2 

where A = ±1. Acting this operator on the wave function (15), one ob- 
tains two subsystems of equations (12) and (13) corresponding to additional 
quantum number A = ±1. The operator (74) obeys the relation = A;^ 
required. It is not difficult to check that all introduced projection operators 
(70), (73) and (74) commute with each other and with the operator of the 
wave equation (20). So, the product of these operators 

Aa'S'(±i),(o)'S'(^i),(o)> (75) 

extracts the states with definite spin, spin projection and additional quantum 
number A = ±1 which is connected with doubling of spin states. It should 
be noted that these states are not states with the definite parity because the 
operator A;s^ mixes states with the opposite parity. The parity operator is 
given by P = r7pr4r4, and it does not commute with the operator A^. To 
have the states with definite energy, we consider Eq. (20) in the momentum 
space: 

eip + mPi + -^Po *(p) = 0, (76) 



where p = PuT^, e — 1 corresponds to positive energy, and £ = —1 - to 
negative energy. It easy to verify, using Eqs. (22), (23), that the projection 
operator 

i (p^^^ i (p^^^ +p^^^) - em 



2m^ 
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ipW UpW — smj ip^^'^ up*^^) — em) 

where p^^^ = Pvl3^\ P^^^ = Pul3i^\ obeys the equahty [M^y = , is a 
solution of Eq. (76), and corresponds to states with spin one (see [12], [ISl)- 
Using Eqs. (21), (56), one obtains the equahties 

f3l'^T, = T,(3^J\ f3^?T,=T,(3l'\ (78) 

It is easy to see with the help of Eqs. (78) that operator (77) commutes with 
the operators (74). Using the properties of the entire algebra (16), and Eqs. 
(21), (70), (73), one may verify the relations [12], [T^ : 

= = 0, = Sf,/^^ = 0, ^fo)P« = Sf,/'^^ = 0, 

p(i)p(o) ^ p(o)p(i) ^ p(T)p(o) = = 0, 

. . . . ' (79) 

p(l)p(l) ^ = p(.0)p{0) ^ p(0)p(0) ^ 

S^o)P^'^ = P^'\ %p(°)=p(°\ 5(V«=J>(^\ 5(V^^=p(^). 

It is checked with the help of Eqs. (79) that the operators (77) commute with 
the operators (75). As a result, from Eqs. (75), (77), we obtain projection 
operators, in the form of matrix-dyads, extracting pure states with spin one, 
spin projection ±1, 0, definite energy, and quantum number A = ±1: 

A« = mW^^i^^dAa = • 

(80) 

The matrix-dyads A^'^\ Aq^-*, correspond to the states with spin one, spin pro- 
jections ±1 and 0, correspondingly. The construction of matrix-dyads for the 
states with spin zero requires finding other solutions of Eq. (76) correspond- 
ing to spin zero, and here it is not considered. The dyad representation (80) 
can be applied to quantum electrodynamic calculations of possesses with the 
presence of fields under consideration. The method of computing the traces 
of 16- dimensional matrix products was developed in [T2|, [IH]- 



18 



8 Canonical quantization of fields 

To apply the canonical quantization [23], we consider the fields "^aIx) in the 
Lagrangian (54) as "coordinates" and the variables do'^Aix) as "velocities". 
The momenta are defined as 



It should be noted that there are no constraints here. With the help of the 
Poisson brackets between "coordinates" "i^Aix) and momenta 7r^(a;), from 
Eq. (81), one obtains 

{^a(x, t),t (^+(x', t)T,)j = 6abS (x - x') . (82) 

— 2 

Using the equality r4 = 1, and the transition to the quantum theory by the 
substitution 

{^,7r} = -i [^,7r], 
where [^,71] = ^E'tt — 7r\E', we arrive at the quantum commutators 

^Aix),ni^')]^^^, = (r4)^^5(x-x'), (83) 

where A, B = 1,2,. ..,16. The quantization is performed here in the same 
manner as in the case m = mo (a=l) [12], [12] • Using Eqs. (15), (19), one 
finds the commutation relations as follows: 

[<^{x), m<^*Q{x%^^, = i(5(x - x'), [<^{x),m<^l{x%^^, = -i6{x - x), 

[m{pk{x), H*Xx%^^, = iSkm^ix - x'), (84) 

[m(pk{x), E^{x')]^^^, = idkn^y^ - x'), 

where = {^/'^)^mnk'fnk, Em = i'-Pmi- Expressiug the components (Pq{x), 
^q{x), iPfj,u from Eqs. (l)-(3), replacing them into Eqs. (84), and taking into 
consideration the commutators of fields, one arrives at the relations 

[V9(x), dQLp* {x')\^^, = -im5(x - x'), (85) 

[(^(x), 9o<^*(x')]i=i/ = im5(x - x'), (86) 

l 

[ifkix), doiPm{x%^^, = —6kmS{^ - x'), (87) 
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2 

[<pk{x), ^olpl{x')]^^^, = -— 5fc„5(x - x'). (88) 

It follows from Eqs. (85)- (88) that vector and pseudoscalar fields have the 
positive metric and pseudovector and scalar fields lead to negative metric 
because of the additional sign (-) 131. The density of the Hamiltonian is 
defined by the relation TC{x) = TrA{x)do'^A{,x) — C{x). With the help of Eq. 
(81), (15), (48), we find 

U{x) = z^+(x)r49o^(x) 
= m ip*^{x)doEm{x) - E^{x)do(fm{x) + ipo{x)doLp{x) - (f*{x)doLpoix) (89) 

-(fl{x)do(p{x) + <^*(x)(9o<^o(a;) - '^ln{x)doH.m{x) + H^{x)do'^rn{x) 

From the expression for four-current J^{x) = i\E'(x)r^\l'(x), one obtains, with 
the aid of Eqs. (15), (48), the total electric charge 

Q = J d^x'^{x)T4^{x) = J d^x^+{x)T4^{x) 
= -im J (fx (p*^{x)E^{x) - E^{x)(p^{x) + (p*Q{x)(p{x) - (p*{x)(po{x) (90) 

-'P*0{X)V{X) + <^*{x)<^o{x) - ^*^{x)Hm{x) + H*^{x)^m{x) . 

The quantizing procedure here is similar to one in quantum electrodynamics 
(QED) for the bispinor fields. But the difference is in the statistics: bispinor 
fields obey the Fermi-Dirac statistics, and 16-component fields \E'(x) under 
consideration, are bosonic fields satisfying the Bose-Einstein statistics. As a 
result, we have here commutators (83) instead of anticommutators in QED. 

It follows from Eqs. (89), (90) that neither energy nor charge of the 
classical fields under consideration have positive values. Therefore in quan- 
tized theory we have to introduce indefinite metric The total space of 
states has two subspaces with positive (Hp) and negative (i7„) square norms. 
According to Eqs. (85)-(88) the vector and pseudoscalar states possess a pos- 
itive square norm, and pseudovector and scalar states have a negative square 
norms. The total space of states represents the direct sum of the two sub- 
spaces Hp and 
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9 Electromagnetic interaction of fields 

The "minimal" electromagnetic interaction is introduced by the substitution 
— > Dy^ = dfj,—ieAfj^, where is the vector-potential of the electromagnetic 
field. The relativistic non-minimal electromagnetic interaction may be taken 
into account by consideration of the operators (see Appendix) : 

where = d^A^ — d^A^ is the strength of the electromagnetic field. Intro- 
ducing phenomenological parameters ki, ^2, i^z in the operators (91), and the 
replacement 5^ D^, one obtains from Eq. (20) the relativistic wave equa- 
tion of the first order for the interaction of bosonic fields under consideration 
with the electromagnetic field: 

2 

T,D, + mPi + ^Po 

m ^ ^ 

1 / - - - ^ 1 ^''^ 

To clear up the physical meaning of constants Ki, K2-, K3 introduced, we 
rewrite Eq. (92) in tensor form. With the help of Eqs. (97)-(99) (see Ap- 
pendix), (15), (19), from Eq. (20), one can obtain the system of equations 
for interacting scalar, vector, pseudoscalar and pseudovector fields with elec- 
tromagnetic field: 

1 

D^,^^Ju. H -i/jq + - (^2 - «i - ^3) ^^i.V'IH = 0, 

2 

jYt Z 

D^,^^, H -^0 + 7; ('«2 - «i + ^3) ^nutp[nu] = 0, 

m 2 

Diy^Pll^u] + D/^lpQ + mV^^ + {k2 + Ki) J^^ui'u + K2- Us) ^tiu^u = 0, 

(93) 

iD^ipl^^] + D^Vo + rnip^ + {k2 + Ki) + i (ki - + ^^3) ^^ui>u = 0, 

D^il^^ - Dnip^ + e^,,apDfiija + rnipif,^] + {k2 + ki) {J^upi^t^p - ^npi^up) 

+ {ki - K2- K3) T/^ulpo + i{Ki- K2 + K3) = 0, 

where dual tensors arc defined by JT^^ = (l/2)£^^a/3J^a/3, V'M = {l/'i)e^uapi^[a(i] 
{^ij.vai3 — {—i)&ixvai3)- Replacing the variables in Eq. (93) with the help of 
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equalities (see Eqs. (l)-(4)) = "V^, V'm = fnf^,, ip[f,u] = <ff,u, '^^i = irriLp^, 
ipQ = —iif, one can express the scalar and pseudoscalar fields from Eqs. (93) 
as follows: 



V = —Df,(p^ + (k2 - Ki - K3) J^^„(p^u, 



(94) 



^ = —Df,(p^ + (k2 - Ki + K3) ^t^u^tiu- 



Excluding the scalar and pseudoscalar fields in Eqs. (93) with the aid of Eqs. 
(94) (with the appropriate replacement of fields), one obtains the system of 
interacting vector, pseudovector and tensor fields: 



+m (Pf, + m {Ki + K2) J^^vVv -m{Ki- K2- K3) ^f^u^iy = 0, 
Du^fMU 2 -^M-^aV'" ~ 7^ 2 (^^2 - Ki + K3) [J^au'fau ) 

+ rn^ip^ + m (ki + K2) J^f.u'fu + m (ki - K2 + Kg) ^^i,(fu = 0, (95) 



+- 



m 



[Ki - K2- K3) ^^luDa'^a - (/«1 -1^2 + 1^3) ^fiuDafa 



2ml 

It follows from Eqs. (95) that constants ki, K2, ^3 are connected with 
the anomalous magnetic moment (AMM) and quadrupole electric moment 
(KEM) of a particle [Hj. The terms in Eqs. (95) containing mo are contributed 
from the scalar and pseudoscalar states. At the limiting case mo — > 00 
(masses of a scalar and pseudoscalar states are infinite), one neglects this 
contribution because the corresponding terms approach to zero. At the par- 
ticular case m = mo, ki = ie/2m, ^2 = ^3 = 0, we arrive at the description 
of vector particles with gyromagnetic ratio g = 2 [22], |IS|- So, the usage 
of the representation of the Lorentz group with high dimension (in our case 
dimension is 16) for the description of vector fields allows us to introduce 
more phenomenological constants of electromagnetic interaction. 
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10 Conclusion 



Equations describing fields which may exist in two mass (m, mo) and spin 
(one and zero) states have been considered. Such fields possess the additional 
symmetry due to the doubling of spin states. This symmetry is an analog 
of dual transformations in electrodynamics. At the equal masses of spin one 
and zero fields, we recover the internal symmetry group 5*0(4, 2) (or locally 
isomorphic to the f/(2, 2) group) investigated in ^2), C^]- This symmetry 
allows us to construct gauge theories with the non-compact groups [T^ . 

The matrices of the relativistic wave equation obtained obey the Dirac 
commutation relations. However, the matrix equation contains the additional 
projection operators which are connected with two mass states of scalar and 
vector fields. At the particular case of equal masses of scalar and vector 
states, m = mo, we arrive at the Dirac-like 16-dimensional equation which 
is equivalent to the Dirac-Kahler equation. Such equation is widely used 
for describing quarks (spin 1/2 fields) on the lattice. However, we consider 
here bosonic fields which correspond to the Bose-Einstein statistics. The 
quantization procedure has been carried out similar to QED, but with the 
commutation relations of fields instead of the anticommutators. Although the 
quantization of such fields requires the introduction of the indefinite metric, 
the theory of fields under consideration can be considered as an effective 
theory. 

The density matrices obtained for spin one fields may be applied for 
calculations of probabilities of different quantum processes in a covariant 
manner. The method considered in allows us to make evaluations of 
physical quantities without using the matrices of first-order equations. 

The considered quaternion form of equations is very convenient for the 
study of internal and Lorentz group symmetries. In these particular cases, 
one may get from Eqs. (64), (66) well-known and important Maxwell, Proca, 
etc. equations in the quaternion form. 

We may speculate that the fields under consideration can be used for 
theoretical schemes describing different objects including sub-quark matter. 
This, however, requires to solve the problem of the physical interpretation of 
quantum fields with indefinite metric. 

We have also studied fields describing particles possessing AMM and 
KEM and, therefore, having the internal structure. Possibly, this scheme 
may be applied for the description of composite systems in nuclear and par- 
ticle physics (see j2H])- 
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APPENDIX. Products of 16-Dimensional Matrices 



Consider products of Petiau-Dufiin-Kcmmer matrices. It is implied that 
all matrices act in 16- dimensional space. Using the properties of the entire 
algebra matrices (16), we obtain 

/3(1)/3(1) _ Jan],[au] , r a,a _ y,n r(}-) R'^-) — p <r^'P 

We write out only nonzero elements of matrices in 16-dimensional space. 
With the help of Eqs. (16), (21), (48), one finds the antisymmetric product 
of 16-dimensional Dirac matrices 



1 1 



-r[^r^] = - (r^r^ - r^r^) = ^m.m _ ^mm + ^m,^ _ ^^.m + ^^.^^ _ 

~x 1 - 
( p"^'^ — p^'"^ 1 -I- pO'[H _ p[H.o J p { _ c-[pM,o 



_P /^p'^.p _ pp.^^ _ pO.lM'^] _L Mfi _ tp i p^AM _ p[H,o \ 

^ ^ (99) 

Using Eqs. (47), we obtain the generators of the Lorentz transformations in 
the 16-dimensional space of wave function (20): 

(100) 
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